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Abstract. In this work a first attempt is made to extract the Equation of State (EoS) 
using experimental results of the meson produced in nuclear collisions at AGS, SPS 
and RHIC energies. The data are confronted to simple thermodynamic expectations 
and lattice results. The experimental data indicate a first order phase transition, with 
a mixed phase stretching energy density between ~ 1 and 3.2 GeV / fm 3 . 



1. Introduction 

Quantum chromodynamics (QCD) predicts a phase transition from hadronic matter 
to the quark and gluon state at a critical temperature T ~ 200 MeV pp. There are 
experimental evidences that nuclear collisions at ultrarelativistic energies induce this 
QCD phase transition and create a new state of matter where the properties of the 
matter are governed by quarks and gluons - such a system is called Quark Gluon Plasma 
(QGP). Current lattice QCD calculations indicate that the order of this transition 
depends on the quark masses, as well as on the baryochemical potential (/ie). In 
particular, a first order phase transition is predicted for large baryochemical potentials, 
while at zero and small baryochemical potentials a "cross over" is predicted. When 
studying nuclear collisions at ultra-relativistic energies the baryochemical potential at 
midrapidity decreases with increasing energy. Therefore, lattice QCD indicates that 
the order of the transition changes with energy (/is) of the collision. However, lattice 
calculations for non-zero hb have large technical uncertainties and the best estimate of 
lattice remains at /is=0. 

In a first order phase transition the pressure increases with increasing temperature, 
until the transition temperature T c is reached, then remain constant during the mixed 
phase, and continue to increase after the end of the mixed phase. Related to this picture 
L. Van Hove [2 J suggested to identify the deconfinement transition in high energy proton- 
antiproton collisions, looking at the variation of average transverse momentum ((pt)) of 
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hadrons as a function of the hadron multiplicity at midrapidity (dN/dy) and searching 
for an increase of (px) followed by a plateau-like behavior and again a subsequent 
increase. The (pr) is expected to reflect the thermal freeze-out temperature (Tf) of 
hadrons and a flow component which can be related to the initial pressure, similarly 
the hadronic multiplicity reflect the entropy density of the system. While the purely 
thermal component of (px) can not be related to the initial temperature which remains 
unmeasurable above T c , however the flow component in the inverse slope can reflect the 
plateau of the pressure during mixed phase. 

It has been observed that the {mj) of pions, kaons and protons as a function of dN c h/dy 
shows a Van-Hove-like behaviour as explained above for a wide range of collision energy 
and for the same centrality [H Hj. Hydrodynamic calculations assuming a first order 
transition could reproduce these data [HJ, [5] . In this work, we study for the first time the 
variation of the inverse slope, T g g extracted from the pt distribution of the <fi -meson as 
a function of the initial energy density, evaluated within the framework of Bjorken's 
hydrodynamical model [7] and y/SNN for energies spanning from AGS, SPS to RHIC. 
The 0-meson is of special interest [6] because of its small hadronic rescattering cross 
section of cr(0N)=lO mb causing it to decouple earlier than other hadrons. Furthermore, 
due to its life time of ~ 45 fm/c, its main decay product (K + K~) suffer less rescattering. 
Experimental results from Au+Au collisions at RHIC energies indicate that the has 
a higher thermal freeze out temperature as compared to pions, kaons and protons. 
In particular their thermal freeze out temperature is within errors compatible with 
the chemical freeze out temperature of hadrons and the critical temperature. Another 
important experimental fact is the observation of a scaling of the elliptic flow of hadrons 
as a function of the transverse momentum when divided with the number of valence 
quarks. This observation is interpreted as an indication that elliptic flow builds in the 
partonic phase and can therefore be reflecting the initial conditions e.g. the initial 
pressure. Therefore the and its flow phenomena are particularly interesting probes 
for studying the EoS and the nature of the phase transition. 

Furthermore, we study here for the first time the inverse slope as a function of the 
initial energy density, e#j. This is an important new feature of such studies, because 
it connects the inverse slope with a parameter characterizing the initial state of the 
collision build up after~ 1 fm/c and which reflects at the same time the collision energy, 
the stopping and the impact parameter of the collision. For example, at a given ^/snn, 
different energy densities could be achieved by changing the colliding nuclei species 
or the impact parameter. The esj is a meaningfull parameter also if equilibrium is 
not reached. It can be directly compared to the critical energy density e c obtained in 
lattice QCD calculations- e c ~ 1 GeV/fm 3 . In the following the data will be analysed 
and confronted to simple thermodynamic expectations which relate to the Van Hove 
signature and to lattice QCD predictions. 
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Figure 1. Tf ff (left ) and (right) as a function of ^/sWn from AGS-SPS to RHIC. 



2. as a probe of the order of the phase transition 

To study the dependence of the inverse slope of the </> meson on the collision energy and 
6b j, we have compiled the inverse slope of the transverse mass spectra (my = ^Jm^ + Pj) 

called "effective temperature" T^j using data at mid-rapidity in the low px domain from 
AGS [8J-SPS P EH m E21 H3f to RHIC PI]. The inverse slope has been estimated 
using an exponential function of the type: = A. exp (yjj)- The use of slightly 

different fit functions leads to differences of the order of 10 MeV, which are small as 
compared to the experimental errors on the inverse slopes [10J. The inverse slope of the 
meson extracted through the above fit includes a thermal component and a non-thermal 
component due to the collective transverse flow [IE]. At low p? (in non-relativistic 
domain, px <C m), T e // = T t h + \m (/?t) 2 , where /3t is the collective transverse velocity 
We have used the above equation to estimate from T^g, as a first approximation. 
The values of {(3t) as a function of ^snn has been obtained from Ref.fTB] and as a 
function of centrality for ^s^n = 200 GeV AuAu collisions from Ref . [T7] . 

In figured] the effective temperature of 0, T^ (left) and the thermal component of 
the inverse slope (right) are shown as a function of the collision energy y / J/v r /v. It is 
observed that from AGS to SPS energies these observables remain almost unchanged, 
showing a plateau-like structure. Going from SPS to RHIC energies, the inverse slope 
eff exn i w ts a sudden jump while an increase is still observed in the component. 
This may be due to an imperfect transverse flow component subtraction at RHIC or 
other effect. It is observed that the is reaching at RHIC values compatible within 
errors with T c . The observed plateau of T^ is a signature of a coexisting phase of 
quarks, gluons and hadrons for a first order phase transition, during which the initial 
pressure remains constant. The subsequent increase of T^ with ^s^n at top RHIC 
energies indicates the end of the mixed phase and the entering into a pure QGP phase. 

Now we study the dependence of the inverse slope and its thermal component on 
the initial energy density, e B j, estimated as e B j = (^f) = (^7/ ( m r) 7^2 where 
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Figure 2. Tf ff and as a function of e Bj .r from AGS-SPS to RHIC. 
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with (m T ), the observable 6b j could be estimated for all centralities and center of mass 
energies at mid-rapidity. 

In figure [2] T^g (left) and (right) are shown as a function of e^.r. Note that 
the formation time, r is model-dependent and in general, in subsequent discussions 
we also assume r ~ 1 fm/c. The above result reflects the properties of equation of 
state. We observe a plateau in T^ stretching between esj ~ 1 and 3.2 GeV//m 3 , and 
increasing suddently above 3.2 GeV//m 3 . This behaviour as already discussed, suggests 
a first order transition, however from figure 2 we can now infer that the mixed phase is 
stretching between 1 and 3.2 GeV/fm 3 . The use of the esj scale, allows us to establish 
here for the first time the €bj range of the mixed phase. More data on the <fi at SBj 
below 1 GeV/ fm 3 are needed to establish the increase of T^ up to 1 GeV//m 3 , as 
seen in other hadrons. The increase of the inverse slope at RHIC energies again indicate 
a pure phase of QGP, as is expected from a first-order phase transition. The thermal 
component shows also a plateau while a smaller increase is still observed at 3.2 
GeV//m 3 . 

The flow component of the inverse slope of the reflecting the initial pressure, it is of 
interest to look directly observables linked to this initial pressure like the transverse flow 
velocity, (3t for all hadrons and the elliptic flow as a function of collision energy. The 
transverse flow velocity (3? shows exactly the same characteristic behaviour as the T^ 
as a function of collision energy, namely increase, a plateau and subsequent increase at 
RHIC [16] . A similar pattern is suggested for the elliptic flow v-i as a function of collision 
energy [18J. 



3. Comparison of <f> data with lattice predictions 

In the following we compare experimental data with lattice QCD predictions. In figure 
[3] the lattice prediction for the energy density is shown as a function of the temperature 
(left). This is compared to esj as a function of Tf h (middle) and T^, (right) extracted 
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Figure 3. Energy density as a function of temperature: lattice plot from Ref. [T] 
(left). Middle and right: Energy density as a function of T t h (middle) and T e ff (right) 
for the c/>. 

from the slope of <f> spectra. First we discuss the middle plot. The lattice prediction is 
for zero baryochemical potential therefore correspondig to the cross over region. The 
figure in the middle (data) is at non-zero \xb and shows the variation of the initial energy 
density with Tf h temperature, which is measured at a later time than the energy density, 
namely at the thermal freeze out of the (p. That is, the two variables in the data are 
a measure of the system at different times, while the lattice estimate is independent of 
time. This temperature is expected to be always below T c and does not reflect the initial 
T, which may exceed T c depending on the collision energy. This plot can be directly 
compared to figure 11 of [15] . where the energy density from data has been studied as 
a function of T, while both x and y axis were estimated at the same time, namely at 
the hadronic chemical freeze out time and at \xb — 0. It is seen that the energy density 
increases approaching the T c from below, from which both the T c and critical exponents 
can be extracted [T9] . 

If the transition occurs at T c ~ 200 MeV as expected, the temperature of <p at the 
thermal freeze out should saturate below and near T c , for all values of the initial density 
(up to infinity), as seen in figure [3] (middle panel). We observe a saturation in the value 
of Tf h near Tc, because Tf h < T c as mentioned earlier. To measure temperatures above 
T c photons and dileptons are useful probes. 

To avoid this saturation and to study an equivalent of the initial temperature, we use in 
the following figure (right) the inverse slope of <ft instead of its temperature at hadronic 
freeze out. Doing that, we take advantage of the possible relationship between collective 
transverse flow and initial pressure. 

We now discuss figure |3] (right panel), where the initial Bjorken energy density is shown 
as a function of the effective temperature of the <fi for collisions from AGS, SPS to RHIC. 
Now both the x and y axis are reflecting parameters at initial times therefore this plot is 
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more appropriate to be compared to the lattice results. The effective temperature of the 
<j) here is a sum of a thermal freeze out temperature, which is expected to be always below 
T c , and a non-thermal component due to transverse flow, which relates to the initial 
pressure and reflects the initial conditions above T c . Therefore the variables in the two 
plots compared here are not exactly the same but they are correlated. Further analysis 
is needed to compare the data to lattice e.g. using exactly the same variables in both 
estimates. A study can be also done involving the elliptic flow v 2 and the transverse flow 
velocity (3t as a function of the initial energy density. Also a hydrodynamic calculation 
of the discussed variables is of interest. The above items are work in progress. 

4. Summary 

In summary, a first analysis of experimental data on the inverse slope parameter of 
mesons as a function of ^snn & energy density from AGS, SPS to RHIC, and their 
comparison to simple thermodynamic expectations, suggest a first order phase transition 
as predicted by lattice QCD at non-zero baryochemical potentials. The mixed phase of 
quarks, gluons and hadrons, is found to stretch between €bj of ~ 1 and 3.2 GeV/fm 3 . 
The latent heat density in a first order phase transition is 6q(T c ) — e#(T c ) = 4B, 
where eg(e#) is the energy density of QGP (hadrons) at T c and B is bag constant. 
With eg(T c ) ~ 3.2 GeV/fm 3 and e#(T c ) ~ 1 GeV/fm 3 , we get a reasonable value for 
B 1 / 4 ~ 250 MeV. If we ignore B, then eg/ 'e# ~ gg/ 'gjj ~ 3.2 also a reasonable number 
comparable to results from lattice QCD [20J, here gQ^gu) is statistical degeneracy of 
QGP (hadrons). The plateau and subsequent increase of the inverse slope parameter 
of above e#j 3.2 GeV/fm 3 is in agreement with data from pions, kaons and protons, 
and is as well observed in the transverse flow velocity (3t reflecting the behaviour of the 
initial pressure. A first attempt to compare data to lattice QCD predictions is made. 

The above results, while supporting the order of the transition predicted by lattice at 
large \lb up to e^j ~ 3.2 GeV/fm 3 , they do not exclude a change of order of the transition 
at smaller baryochemical pontials. Which points to further work, towards mapping 
out the order of the QCD phase transition as a function of energy and baryochemical 
potential. 
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